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Abstract 

We have measured the thermal conductivity of the heavy-fermion superconductor CeCoIns in 
the vicinity of the upper critical field, with the magnetic field perpendicular to the c axis. Thermal 
conductivity displays a discontinuous jump at the superconducting phase boundary below critical 
temperature Tq ~ 1 K, indicating a change from a second to first order transition and confirming 
the recent results of specific heat measurements on CeCoIns. In addition, the thermal conductivity 
data as a function of field display a kink at a field below the superconducting critical field, 
which closely coincides with the recently discovered anomaly in specific heat, tentatively identified 
with the appearance of the spatially inhomogeneous Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) 
superconducting state. Our results indicate that the thermal conductivity is enhanced within the 
FFLO state, and call for further theoretical investigations of the order parameter's real space 
structure (and, in particular, the structure of vortices) and of the thermal transport within the 
inhomogeneous FFLO state. 

PACS numbers: 74.70.Tx, 7L27.+a, 74.25.Fy, 75.40.Cx 
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I. INTRODUCTION 



as well as the stable spatial 
ig . are subjects of intense 



Over the last several years there has been renewed interest in the spatially inhomogeneuos 
Fulde-Ferrell-Larkin-Ovchinnikov(FFLO) state. The FFLO state was predicted as early as 
the mid-1960s to occur in a clean Type II superconductor in high magnetic fields, 

when the Zeeman energy becomes comparable to the condensation energy. Then, Pauli 
limiting Q plays an important role in defining both the superconducting critical field ifc2 
and the temperature Tq below which the FFLO state is expected to appear Within the 
FFLO state spin up and spin down electrons of a spin-singlet superconductor can only stay 
bound if the Cooper pair has a finite momentum. As a result, the FFLO state is formed with 
a spatially oscillating order parameter. The exact description of the corresponding phase 
diagram for both s- and d-wave superconductors, js, 6, 7, 
structures in 2D and 3D in the presence of vortices 5|, |9 , 
theoretical investigations. 

In spite of the straightforward nature of the theoretical prediction, the experimental 
observation of the FFLO state has turned out to be a difficult task. In fact, very few 
superconductors fulfil the necessary conditions for the formation of an FFLO state. The 
relative importance of the Pauli and orbital limiting can be described by the so-called Maki 
parameter a = v2-^. if^g is the orbital limiting field due to the kinetic energy of the 
superconducting currents around the vortex cores, commonly derived from the slop e of the 
experimentally determined H — T phase boundary at Tc as H^2 = 0-'''^^lrc 'l-^- Hp = 
^/2Ao / gfiB is the Pauli limiting field due to the potential energy of electrons' spins (Zeeman 
energy). Here Aq is the zero temperature value of the superconducting gap, g is the electron's 
effective g-factor, and fiB is the electron's Bohr magneton [^. Within the calculation of 
Ref. Q, 

a must be greater than 1.8 for the FFLO state to be realized. 
There are several classes of materials that are traditionally thought of as potential candi- 
dates for the formation of the FFLO states. These include low dimensional organic supercon- 
ductors and heavy fermion superconductors. The low dimensional organic superconductors 
are promising because when the field is applied within the conducting planes of a 2D super- 
conductor, the orbital limiting is suppressed entirely, as the diamagnetic screening currents 
can only fiow within the plane. In such case Pauli limiting determines the critical field 
Hc2 = Hp, the Maki parameter a = oo, and the FFLO state should be stabilized below the 
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critical temperature Tq ~ 0.55Tc for magnetic field close to Hc2- This straightforward pre- 
diction led to a number of experimental investigations of the superconducting properties of 
lower-dimensional organic superconductors. Several investigators suggested the existence of 
FFLO =tate, e.g. b^ed on the superconducting phase d.ag.an, Q o. the n,agneto-thcnnal 
transport properties |14| . 

Heavy fermion superconductors are also attractive because of their potentially large values 
of the Maki parameter. Here the heavy electron masses lead to low Fermi velocities of the 
quasiparticles and in turn to relatively ineffective orbital limiting, or large H^2- For this 
reason heavy fermion materials have had their deserved share of attention, and the features 

n 

in the magnetization of CeRu2 [15[ and the phase diagram of UPd2Al3 |16|] were taken as 
possible signatures of the FFLO states. 

However, and not for the lack of effort, there is to date no accepted definitive proof of 
the existence of the FFLO state in any of the systems described above. In this paper we 
describe the magneto-thermal studies of a new and very strong candidate to possess the 



FFLO state, the heavy-fermion superconductor CeCoIuR. 

CeCoIns is a clean d-wave superconductor, with a Tc of 2.3K, the highest 



among the Ce-based heavy fermions. It exhibits a layered structure of alternating Celus and 
CoIn2 planes, suggesting a possible quasi-2D electronic nature of this compound. This is 
supported by the experimental observation of a quasi-cylindrical sheet in the Fermi surface 
of CeCoIns via de Haas- van Alphen studies j^l • The estimated Maki parameter a is about 
3.6 for the field perpendicular to the Celns planes {H \\ c) j2]|, and close to 4.5 for the field- 
in-the-plane orientation {H ± c) 22]. Thus, CeCoIns is a good candidate for the formation 



of an FFLO state. As new results accumulate, there is a growing evidence that the FFLO 
state may indeed be realized in CeCoIns. First, the superconducting phase transition at low 
temperatures becomes first order, which is manifested by the sharp specific heat anomaly 



at Tc in the specific heat for both H 



steps in magnetostriction J21|, |2J] and in magnetization 



and H ± c m 



orientations. In addition, 
and a step in the thermal 



conductivity for H \\ c are observed. The change of the superconducting anomaly 
from second to first order j2l| was interpreted as a realization of the Maki scenario, which 
attributes this change to a strong Pauli limiting effect in a Type II superconductor j27|, |28[ . 
When the field was applied within the a — b plane {H _L c ), a second anomaly in the specific 
heat was observed within the superconducting state |2^ |2^ , indicating a phase transition 
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into a new superconducting state, tentatively identified as the FFLO state in CeCoIns. 
In addition, steps in magnetization of CeCoIns were observed by Radovan et al. [2^, and 
were interpreted as an indication of the multi-quantum vortices expected under certain 
circumstances for the 2D superconductors within the FFLO state 0, liol . \\\ . The validity 



ImI- On 



30, mi. On the theoretical front, recent 



of such interpretation is at present under debate 
analysis of a linearly increasing Hc2 at the lowest temperatures suggested that this too 
can be accounted for within an FFLO scenario for CeCoIns. 

While these results make the FFLO scenario a very appealing one for CeCoIus, there is 
no clear evidence so far for spatially inhomogeneous superconductivity in the second low 
temperature phase. A recent study revealed an increased penetration depth at the lower 
transition, which was interpreted as a decrease of the superfiuid density due to the formation 
of the FFLO state [32]. An ultrasound investigation of the high field state revealed the 
decrease of the sound velocity from that in the vortex state, which was also presented in 
support of the FFLO nature of that state 

Here we present our results of magneto-thermal transport measurements in CeCoIns with 
the field applied within the Celns planes {H ± c). The LO structure, which emerged from 
the original theoretical work P], is a collections of periodically spaced planes of nodes of 
the superconducting order parameter that are perpendicular to the direction of the applied 
field. The LO order parameter is described as ipif) = t/^q cos(gr'), oscillating in space along 
the direction of vector g* || H a.s illustrated in Fig. In recent years, thermal conductivity 
was used effectively to probe the anisotropy of the order parameter in unconventional super- 
conductors, specifically, the structure of the nodes in fc-space. This is due to the fact that 
normal quasi-particles, which are easily excited along the nodal directions, do carry heat. 



whereas the superconducting background does not [19|, |3J, |35[- The inherent anisotropy 
of the LO state makes thermal conductivity an attractive tool for identifying it. For the 
vast majority of clean superconductors thermal conductivity drops as the sample enters the 
superconducting state due to the opening of the superconducting gap over the entire Fermi 
surface, and the resulting rapid decrease of the number of normal quasiparticles which carry 
heat. One would then expect the thermal conductivity within the nodal planes to be higher 
than in the rest of the sample. Thermal conductivity would then be larger when the heat 
current Q flows along the nodal planes, and perpendicular to the applied magnetic field 
(Q _L H), than when the heat flow is parallel to the magnetic field {Q \\ H). Note that this 
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anisotropy is of the opposite sign from that due to the vortices, with higher thermal conduc- 
tivity along the vortices {Q \\ H), |36| and therefore the two contributions should be easily 
differentiated, especially if the contribution due to the 2D planes turns out to dominate that 
from the ID vortices. This picture unfortunately turns out to be too simplistic for the case 
of CeCoIus and is complicated by several effects described below. 

Motivated by the idea described above, we measured thermal conductivity in CeCoIns at 
low temperatures, in the vicinity of the upper critical field, with the magnetic field oriented 
in-plane, using a dilution refrigerator in the 20 T magnet of the NHMFL facility at LANL. 
The sample, a needle-like single crystal with dimensions of(2.18 mm, 0.28 mm, 0.064 mm), 
was flux-grown at LANL as described in Ref. After a chemical etch and polishing to 
remove the residual free indium, the sample had resistivity of 3.2 fiQ — cm at 4.2 K and a 
RRR of p(4°2^j = 9.4. The experimental set-up for thermal conductivity consists of a heater 
attached to one end of the sample, and two Ru02 thermometers in thermal contact with the 
sample at two points along its length. A DC heat current flows along the [100] direction of 
the sample (along its longest dimension). The resulting temperatures of both thermometers 
are measured using an LR700 resistance bridge. The measurements were performed with 
magnetic field applied parallel and perpendicular {H \\ [010]) to the heat current. The 
two thermometers were calibrated at each field against a reference thermometer placed in a 
field-free region. 



II. THERMAL CONDUCTIVITY IN THE NORMAL STATE AND ZERO FIELD 

Fig. 121 shows the thermal conductivity of CeCoIns as a function of temperature up to 
2.5 K at zero field and at 12.5 T, for field oriented parallel and perpendicular to the heat 
current. One notices a substantial drop in thermal conductivity induced by the applied 
field above Tc as well gnificant difference between the two field orientations. It 

is possible to account for the normal state anisotropy by considering magnetoresistance of 
CeCoIus, displayed in the inset of Fig. El The longitudinal and transverse magnetoresistances 
at 2.5K and 9T, with an in-plane current, are 37% and 60% respectively. This gives a 
difference in magnetoresistance {p± — p\\)/ p(0) of 23% between the field oriented parallel and 
perpendicular to the electrical current. If we determine the difference of thermal conductivity 
between both field orientations, {k\\ — k,±)/k,{0) at 2.5 K and 12.5 T we find a value of 21%, 
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very close to the anisotropy of magnetoresistance. Thus, the anisotropy of the heat transport 
in the normal state can be accounted for by the difference in the quasiparticle scattering 
for the two field orientations. This is not surprising because the quasiparticle contribution 
dominates the heat transport in CeCoIus Therefore, in order to highlight the differences 
between the superconducting states for the two orientations of the magnetic field studied, 
in what follows, we often present the thermal conductivity data scaled by the values in the 
normal state at 12.5 T. 

The normal state thermal conductivity for the field of 12.5 T, just above the supercon- 
ducting critical field of 12 T, is displayed in Fig. |21^d) as |; vs. T on a log-log plot. ^ 
appears to diverge, reflecting the possible presence of the quantum critical point (QCP) in 
CeCoIns, suggested by both specific heat and resistivity measurements j^. Similar behav- 



ior with t 



le QCP lying very close to the superconducting critical field was observed for H 



[001] |38 

It is interesting to compare the zero-field thermal conductivity measurements on the 
present sample to the previously published data in particular the zero-temperature 
limit of ^. The zero-field data are displayed in Fig. Efc) as |; vs. T^, and appear to 
extrapolate to the value of a few tenths of W/Km, close to the value obtained for the original 
sample. This agreement supports the existence of the universal, independent of the impurity 
concentration, limit of ^ in CeCoIus, and the original interpretation of the low temperature 
thermal transport reflecting the presence of lines of nodes in the superconducting energy 
gap of CeCoIns Q. 

III. THERMAL CONDUCTIVITY IN THE VORTEX STATE 

The thermal conductivity data in the low temperature - high field part of the phase dia- 
gram are displayed in Fig. IHl The transition to the normal state is marked by a pronounced 
jump in the thermal conductivity at the lowest temperatures. The jump in thermal conduc- 
tivity confirms the first order nature of the superconducting transition, reported previously 
on the basis of the specific heat measurements js^. The first order nature of the super- 
conducting transition for H \\ [00 1] was deduced on the basis of the thermal conductivity 
measurements by Izawa et al. jl9| . 

The normalized thermal conductivity for both parallel and perpendicular field orienta- 
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tions is depicted in Fig.|H^a,b) as a function of temperature for several values of field between 
50 mK and 2.5 K. The inset in Fig. |3] shows the data after subtraction of the normal state 
12.5 T data. The enhancement of thermal conductivity below T^ at zero field due to an 
increased quasiparticle mean free path, shown in Fig. |2^a), can still be clearly resolved at 9 
T. An apparent small rise in thermal conductivity at the superconducting transition at 10 T 
and 10.4 T is due to the competition between the increase of the thermal conductivity in the 
normal state with increasing magnetic field and a drop in thermal conductivity as the system 
becomes superconducting when magnetic field is swept (see Fig. EI). When temperature is 
swept, the effect of the increasing k in the normal state prevails at 10 T and 10.4 T. The 
data for 10.8 T show a pronounced step (similar to the data for the field sweeps) associated 
with the first order nature of the superconducting transition, which dominates at this field. 

Fig. |S1 displays the thermal conductivity normalized by the value in the normal state at 
Tc at several temperatures for the two field orientations between 8 T and 12 T. The thermal 
conductivity decreases with increasing field in the mixed state and increases slightly in 
the normal state. The absolute slope of the thermal conductivity vs. magnetic field in 
the vortex state is reduced as the temperature is lowered. This is likely a result of the 
competition between the effects of temperature and magnetic field. Magnetic field has 
two competing effects on the quasiparticle heat transport. First is the Volovik effect, or 
doppler shift of the quasiparticle energies, which results in the finite density of states in 



d-wave superconductors 



42]. This leads to the increase in thermal conductivity with 



the increasing magnetic field. The second effect is due to the quasiparticle scattering from 
vortices. As the field is increased, so is the number of scattering vortices, and this has the 
effect of decreasing thermal conductivity in higher magnetic field. At high temperature the 
number of quasiparticles is largely determined by temperature, and the contribution from 
Volovik effect loses its significance. The vortex scattering effect then dominates thermal 
transport, resulting in the decrease of thermal conductivity with magnetic field. On the 
other hand, at low temperature the Volovik effect dominates the thermal broadening, and 
efficiently competes with the reduction of thermal conductivity due to the vortex scattering. 
In CeCoIus this results in a slower decrease of thermal conductivity with increasing magnetic 
field at lower temperatures, as displayed in Fig. We do not reach the \/H regime observed 
in YBCO [4^ or BSCCO [4^] high temperature superconductors. Perhaps CeCoIns must be 
cooled to lower temperature to observe this field dependence due to its very small impurity 
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band width In the temperature and field range studied, the quasiparticle scattering in 
CeCoIns appears to be dominated by vortices, rather than the impurities. 

The overall values of thermal conductivity ioi H || J orientation are higher than that 
for H J- J. At 0.21 K and 0.27 K the difference is close to 12%, and at higher tempera- 
tures it is significantly reduced (only 4% at 0.57 K and 2% at 0.81 K). This anisotropy of 
thermal conductivity is due to the vortex scattering of the quasiparticles. More precisely, 
in a semi-classical approach the scattering off the vortex is maximal when the quasiparticle 
velocity is perpendicular to it, resulting in a lower thermal conductivity when the field is 
perpendicular to the heat current. This anisotropy is naturally expected to vanish at Hc2 



when the vortices overlap |36|. The same calculation predicts that the anisotropy of k will 
change sign at the T << Tc limit, where the excitation of the quasiparticles perpendicular 
to the vortices (Volovik effect) becomes the dominant effect and enhances the heat current 
in the direction perpendicular to the magnetic field. Our data show that the anisotropy of 
thermal conductivity in CeCoIns is growing to the lowest temperature measured, once again 
indicating the need to go to lower temperature to test the theoretical prediction. To our 
knowledge there are no similar calculations for the d-wave case. More theoretical work is 
needed to help us understand the magnetothermal transport in the vortex state of CeCoIus. 

IV. THE SUPERCONDUCTING PHASE DIAGRAM OF CECOIN5. 

The upper critical field Hc2, determined from both the temperature and the field sweeps, 
is displayed in Fig. [Altogether with the critical field deduced from the specific heat. Note 
that there is no difference in Hc2 for the field parallel and perpendicular to the heat current, 
as expected for a tetragonal compound, since the additional in-plane anisotropy due to the 
d-wave gap is zero upon the 90 deg rotation. The shape of Hc2 as well as the first order 
nature of the transition at low temperatures is in agreement with previous reports 0, |3| ■ 
The strong temperature dependence of B.c2 at low temperatures, as opposed to the saturated 
behavior expected in the BCS theory, is an important observation supporting the existence 
of the FFLO state in CeCoInc: 1221. In fact, this is a common feature for a number of other 

n 

superconductors suggested to be in the Pauli-paramagnetic limit, such as UBeis and 
K - {BEDT - TTF)2Cu{CSN)2 H- 

The amplitude of the jump is comparable for both field orientations and decreases almost 
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linearly with increasing temperature, as shown in the inset of fig. El This is in contrast to the 
sharp decrease in the temperature step, associated with the first order transition, observed 
at the critical point in magnetocaloric measurements for the field along the c-axis j2l|. It 
would be difficult to locate precisely the critical point where the order of the superconduct- 
ing phase transition changes from second to first, for H || J, because the transition occurs 
gradually. Nevertheless, the critical point is consistent with that determined from specific 
heat measurements 1231. 



V. FULDE-FERRELL-LARKIN-OVCHINNIKOV TRANSITION. 

An additional feature in the k vs. H curves displayed in Fig. El is resolved at the low- 
est temperatures. A kink appears in the data at a field Hk for H || J, and the thermal 
conductivity is nearly constant between and ifc2 for temperature below 0.27 K. Given 
the sharpness of the jump at Hc2 and the first order nature of the transition, this feature is 
clearly distinct from the Hc2 anomaly. This anomaly is not present in the data at 0.57 K 
and above, where the thermal conductivity continues to decrease up to Hc2. We therefore 
interpret the region of the fiat i^{H) preceding the sharp jump at Hc2 as an enhancement of 
thermal conductivity. Further, decreases as the temperature is reduced. Fig. El displays 
iffc in the H — T plane, together with the phase diagram deduced from the specific heat 
measurements The points coincide well with the second low temperature phase 

transition line found in the specific heat below Hc2. The absence of such a fiat portion at 
high temperatures, where the thermal conductivity is monotonously decreasing up to Hc2, 
and the good agreement of Hjt with the specific heat anomaly lead us to speculate that the 
enhancement of thermal conductivity is due to the formation of the second low temperature 
superconducting state in CeCoIus, identified as a potential FFLO state 23, |2^. Features 



in the thermal conductivity of an organic superconductor near the superconducting critical 



field were also interpreted as possible signatures of the FFLO state 

The enhancement of thermal conductivity between and ifc2 manifests itself clearly 
in the H || J data. It is more difficult to make a definitive statement for the data in the 
H J- J geometry since data below Hc2 is rounded and a gradual rise is present in the data 
for higher temperature, outside of the FFLO phase. 

The enhancement of thermal conductivity for the H \\ J orientation is contrary to our 
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original simple-minded expectation, since in this geometry the nodal planes in the LO phase 
are perpendicular to the direction of the heat current, and therefore would not be expected 
to enhance thermal conductivity. There are several possible explanations of our results. 
Recent theoretical work ^ suggests that the lowest energy state is not a pure LO state 
with a single modulation wavevector Q, but a modified LO state with a combination of 
three modulation wavevectors. If so, one would not expect an additional anisotropy with 
respect to the direction of the magnetic field due to FFLO nodal planes. The contribution 
to thermal conductivity from the nodal planes of the LO state is not a priori dominant over 
that from the vortices, and must be investigated theoretically [4^. Another scenario 
suggests that the bottle neck for the heat transport along the field direction is the vortex 
cores. One would expect the structure of the vortices to be modulated by the nodal planes. 
The vortex cores' size might increase at the nodal planes, reducing the bottle neck and 

ay between the vortex and 



The resulting spatial 



leading to the enhancement of thermal conductivity. The interp 
FFLO state was theoretically studied for 2D superconductors 
structures can be alternating nodal planes and lines of vortices, or more intricate structures, 
depending on the Landau level quantization number of the order parameter. More theoretical 
work on the vortex structure within the FFLO state for 3D superconductors, and its effect 
on the thermal transport in particular, is called for and should clarify whether such new 
structures can account for the observed enhancement of the thermal conductivity in the low 
temperature superconducting state. 



VI. CONCLUSION 



Thermal transport is a powerful probe of a superconducting state. We performed thermal 
conductivity measurements to investigate the properties of both the vortex state and the 
second low temperature phase of CeCoIns, with the magnetic field applied within the a — b 
plane of this tetragonal compound. Our data demonstrate that the superconducting phase 
transition becomes first order between roughly 10 T and the superconducting critical field, 
in accord with previous specific heat measurements, indicating the importance of the Pauli 
limiting effect in CeCoIns. In addition, we observed a kink in thermal conductivity for the 
field parallel to the direction of the heat current, coincident with the phase transition in the 
second low temperature state of CeCoIns, suggested previously to be an FFLO state [3, 
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29|. Thermal transport within the FFLO state at present remains unexplored theoretically, 
and the observed enhancement of thermal conductivity within the FFLO state of CeCoIns 
is puzzling. Our experimental results present a new challenge for the understanding of 
inhomogeneous superconductivity. 

VII. ACKNOWLEDGEMENTS 

The authors are grateful to I.Vekhter, L.Bulaevskii, P.Fulde, F.Ronning, M.Tanatar for 
fruitful discussions. 



* Present address: Hochfeldlabor Dresden Forschungszentrum Rossendorf Postfach 51 01 19 
01314, Dresden Germany. 
[1] P. Fulde and R. A. Ferrell, Physical Review 135, A550 (1964). 

[2] A. I. Larkin and Y. N. Ovchinnikov, J. Exptl. Theoret. Phys. (USSR) 47, 1136 (1964), [Sov. 

Phys. JETP 20, 762, (1965).]. 
[3] A. M. Clogston, Phys. Rev. Lett. 2, 9 (1962). 

[4] L. W. Gruenberg and L. Gunther, Phys. Rev. Lett. 16, 996 (1966). 

[5] M. Houzet and A. Buzdin, Phys. Rev. B 63, 184521 (2001). 

[6] D. F. Agterberg and K. Yang, J. Phys. Condens. Matter 13, 9259 (2001). 

[7] M. Tachiki, S. Takahashi, P. Gegenwart, M. Weiden, M. Lang, C. Geibel, F. Steglich, R. Mod- 
ler, C. Paulsen, and Y. Onuki, Z. Phys. B 100, 369 (1996). 

[8] R. Combescot and C. Mora, cond-mat/0311042. 

[9] A. L Buzdin and J. P. Prison, Phys. Lett. A 218, 359 (1996). 
[10] H. Shimahara and D. Rainer, J. Phys. Soc. Japan 66, 3591 (1997). 
[11] U. Klein, D. Rainer, and H. Shimahara, J. Low Temp. Phys. 118, 91 (2000). 
[12] N. R. Werthamer, E. Helfand, and P. C. Hohenberg, Phys. Rev. 147, 295 (1966). 
[13] J. Singleton, J. A. Symington, M. S. Nam, A. Ardavan, M. Kurmoo, and P. Day, J. Phys. 

Cond. Matter 12, L641 (2000). 
[14] M. A. Tanatar, T. Ishiguro, H. Tanaka, and H. Kobayashi, Phys. Rev. B 66, 134503 (2002). 



11 



[15] A. Huxley, C. Paulson, O. Laborde, J. Ttiolence, D. Sanchez, A. Junod, and R. Calemczuk, 

Journal of Physics: Condensed Matter 5, 7709 (1993). 
[16] K. Gloos, R. Modler, H. Schimanski, C. Bredl, C. Geibel, F. Steglich, A. Buzdin, N. Sato, 

and T. Komatsubara, Physical Review Letters 70, 501 (1993). 
[17] C. Petrovic, P. G. Pagliuso, M. F. Hundley, R. Movshovich, J. L. Sarrao, J. D. Thompson, 

and Z. Fisk, J. Phys. Condens. Matter 13, L337 (2001). 
[18] R. Movshovich, M. Jaime, J. D. Thompson, C. Petrovic, Z. Fisk, P. G. Pagliuso, and J. L. 

Sarrao, Phys. Rev. Lett. 86, 5152 (2001). 
[19] K. Izawa, H. Yamaguchi, Y. Matsuda, H. Shishido, R. Settai, and Y. Onuki, Phys. Rev. Lett. 

87, 057002 (2001). 

[20] D. Hall, E. C. Palm, T. P. Murphy, S. W. Tozer, Z. Fisk, U. Alver, R. G. Goodrich, J. L. 

Sarrao, P. G. Pagliuso, and T. Ebihara, Phys. Rev. B 64, 212508 (2001). 
[21] A. Bianchi, R. Movshovich, N. Oeschler, P. Gegenwart, F. Steglich, J. D. Thompson, P. G. 

Pagliuso, and J. L. Sarrao, Phys. Rev. Lett. 89, 137002 (2002). 
[22] H. Won, K. Maki, S. Haas, N. Oeschler, F. Weickert, and P. Gegenwart, cond-mat/0306548. 
[23] A. Bianchi, R. Movshovich, C. Capan, P. G. Pagliuso, and J. L. Sarrao, Phys. Rev. Lett. 91, 

187004 (2003). 

[24] T. Takcuchi, H. Shishido, S. Ikcda, R. Settai, Y. Haga, and Y. Onuki, J. Phys. Cond. Matter 
14, L261 (2002). 

[25] T. Tayama, A. Harita, T. Sakakibara, Y. Haga, H. Shishido, R. Settai, and Y. Onuki, Phys. 

Rev. B 65, 180504 (2002). 
[26] T. P. Murphy, D. Hall, E. C. Palm, S. W. Tozer, C. Petrovic, Z. Fisk, R. G. Goodrich, 

P. Pagliuso, J. L. Sarrao, and J. D. Thompson, Phys. Rev. B 65, 100514 (2002). 
[27] K. Maki and T. Tsuneto, Progress of Theretical Physics 31, 945 (1964). 
[28] K. Maki, Phys. Rev. 148, 362 (1966). 

[29] H. A. Radovan, N. A. Fortune, T. P. Murphy, S. T. Hannahs, E. C. Palm, S. W. Tozer, and 
D. Hall, Nature 425, 51 (2003). 

[30] R. Movshovich, A. Bianchi, C. Capan, M. Jaime, and R. G. Goodrich, Nature, in print. 
[31] H. A. Radovan et ai, Nature, in print. 

[32] C. Martin, C. C. Agosta, S. W. Tozer, H. A. Radovan, E. C. Palm, T. P. Murphy, and J. L. 
Sarrao, cond-mat/0309125. 

12 



[33] T. Watanabe, Y. Kasahara, K. Izawa, T. Sakakibara, C.J. van der Beek, T. Hanaguri, H. 

Shishido, R. Settai, Y. Onuki, and Y. Matsuda, |cond-mat/0312062l 
[34] H. Aubin, K. Behnia, M. Ribault, R. Gagnon, and L. Taillefer, Phys. Rev. Lett. 78, 2624 

(1997). 

[35] I. Vekhter and R J. Hirschfeld, Physica C 341, 1947 (2000). 
[36] K. Maki, Phys. Rev. 158, 397 (1967). 
[37] C. Capan et al., unpublished. 

[38] V. A. Sidorov, M. Nicklas, P. G. Paghuso, J. L. Sarrao, Y. Bang, A. V. Balatsky, and J. D. 

Thompson, Phys. Rev. Lett. 89, 157004 (2002). 
[39] J. Paghone, M. A. Tanatar, D. G. H. ad Etienne Boaknin adn R. W. Hill, F. Ronning, 

M. Sutherland, L. Taillefer, C. Petrovic, and P. C. Canfield, Phys. Rev. Lett. 91, 246405 

(2003). 

[40] A. Bianchi, R. Movshovich, I. Vekhter, P. G. Pagliuso, and J. L. Sarrao, Phys. Rev. Lett. 91, 

257001 (2003). 
[41] G. E. Volovik, JETP Lett. 58, 469 (1993). 
[42] L Vekhter and A. Houghton, Phys. Rev. Lett. 83, 4626 (1999). 

[43] M. Chiao, R. W. Hill, C. Lupien, B. Popic, R. Gagnon, and L. Taillefer, Phys. Rev. Lett. 82, 
2943 (1999). 

[44] H. Aubin, K. Behnia, S. Ooi, and T. Tamegai, Phys. Rev. Lett. 82, 624 (1999). 
[45] H. A. Radovan, R. J. Zieve, J. S. Kim, and G. R. Stewart, |cond-mat /030'645T] 
[46] F. Zuo, J. S. Brooks, and et al., Phys. Rev. B 61, 750 (2000). 
[47] R. Combescot, C. Mora, lco^mat/0311042j 
[48] L Vekhter, private communication. 
[49] L. Boulaevskii, private communication. 



13 



H 



[010] 



[100] 




D 



(r) . h| 




FIG. 1: Illustration of the vortex structure (solid lines) and the FFLO modulation (dashed lines) 
with the field parallel (top) and perpendicular (bottom) to the heat current. 
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T(K) T(K) 

FIG. 2: (a): Thermal conductivity k vs. temperature of CeCoIn5 at zero field and 12.5 T, with the 
field parallel (•) and perpendicular (□) to the heat current. Inset (b): In-plane magnetoresistance 
of CeCoIn5 with field perpendicular (top curve) and parallel (bottom curve) to the current, (c): 
Thermal conductivity of the present sample (A) and the sample used in Ref. [18] (★), normalized 
to the value at Tc of k = 2 W/Km. (d) Thermal conductivity divided by temperature for the field 
of 12 T parallel (•) and perpendicular (□) to the heat current, showing diverging k/T as T — 0. 
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FIG. 3: Thermal conductivity vs. magnetic field of CeCoIns between 8T and 12. 5T. Left Panel: 
H II J. Right Panel: H ^J. (*) 1.07 K, (y) 0.81 K, (•) 0.67 K, (o) 0.57 K, (□) 0.27 K, (o) 0.24 
K, (A) 0.21 K, (*) 0.15 K. 
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FIG. 4: Temperature dependence of thermal conductivity k in CeCoIns between 50 mK and 2.5 
K. (•) H II J; (□) H J- J. Left panel: 10 T; right panel: 10.8 T; inset: k — kn vs. temperature at 
9 T (o), 10 T (v), 10.4 T and 10.8 T (A). The values of k at 12.5 T have been used as kn, 
and the data for different fields have been shifted vertically for clarity. 
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FIG. 5: Normalized thermal conductivity vs. magnetic field of CeCoIns for H \\ J (•) and H J. J 
(□). Upper left: 0.21 K, upper right: 0.27 K, lower left: 0.57 K, lower right: 0.81 K. The data 
for both orientations have been normalized by the corresponding normal state values at 12.5 T. 
Arrows indicate positions of the kinks for H \\ J. 



18 



12- 



11- 



CM 



10- 



9- 



0.0 



0.18 



0.12 



• H //J 

□ H_LJ 



0.2 0.4 0.6 0.8 1.0 1.2 
T(K) 



r 

-A 




★ 

V 



froiT^ ^ 
from C 

H, from K 

k 



0.5 



1.0 



T(K) 



1.5 



FIG. 6: Magnetic field vs temperature phase diagram of CeCoIns deduced from thermal con- 
ductivity and specific heat measurements, (solid line) Hc2 and (*) TpFLO are deduced from the 
specific heat data of Ref [23]. Hc2 (A) and (v) obtained from the thermal conductivity data 
(details are in the text). Arrow indicates the critical temperature where superconducting transi- 
tion changes from second to first order. Inset: the size of the jump in k at the superconducting 
transition for the field sweeps for (•) H \\ J and (□) H ± J. 
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